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The equilibrium structure and functional properties exhibited by brownmillerite oxides, a family 
of perovskite-derived structures with alternating layers of BOq octahedra and BO 4 tetrahedra, viz., 
ordered arrangements of oxygen vacancies, is dependent on a variety of competing crystal-chemistry 
factors. We use electronic structure calculations to disentangle the complex interactions in two 
ferrates, Sr 2 Fe 205 and Ca 2 Fe 205 , relating the stability of the equilibrium (strain-free) and thin film 
structures to both previously identified and newly herein proposed descriptors. We show that cation 
size and intralayer separation of the tetrahedral chains provide key contributions to the preferred 
ground state. We show the bulk ground state structure is retained in the ferrates over a range of 
strain values; however, a change in the orientation of the tetrahedral chains, i.e., a perpendicular 
orientation of the vacancies relative to the substrate, is stabilized in the compressive region. The 
structure stability under strain is largely governed by maximizing the intraplane separation of the 
‘dipoles’ generated from rotations of the Fe 04 tetrahedra. Lastly, we hnd that the electronic band 
gap is strongly influenced by strain, manifesting as an unanticipated asymmetric-vacancy alignment 
dependent response. This atomistic understanding establishes a practical route for the design of 
novel functional electronic materials in thin film geometries. 


I. INTRODUCTION 

The family of brownmillerite oxides (general formula 
ABO2.5 or A2B2O5) are highly studied for use in ionic con¬ 
ducting and anion insertion applications.Their struc¬ 
ture type can be thought of as an ABO^ perovskite with 
one-sixth of the oxygen atoms removed, creating parallel 
rows of ordered anion vacancies along the [ 110 ] crystal¬ 
lographic direction; this results in alternating layers of 
corner-connected BO4 tetrahedra and BOq octahedra 
(Figure la). While cooperative distortions and octahe¬ 
dral rotations are well understood in perovskites, the 
presence of tetrahedral layers in the brownmillerites adds 
additional structural complexity and possible degrees of 
freedom for both structure and electronic function design. 
In addition to rotations of the octahedra, each tetrahe¬ 
dral chain can ‘twist’ in a “left-handed” or “right-handed” 
sense, resulting in two different types of chains related by 
mirror symmetry (Figure lb). Furthermore, these differ¬ 
ent chains can be ordered relative to each other within the 
brownmillerite unit cell, resulting in a variety of structures 
displaying different space group symmetries. 

When the tetrahedra and octahedra are not rotated, 
the structure displays the Imma space group (Figure 2 a); 
this aristotype may be used as a high symmetry reference 
phase for which subsequent structural analyses are made. 
Alternatively, the Imma structure can be imagined as 
having disordered left- and right-handed chains {i.e. the 
chains display either incommensurate ordering or no long 
range order); some structures, such as Sr2CoFe05 and 
Sr2MnGe05,^’^^ display this phase at ambient conditions, 
while others become disordered at high temperature (in- 



FIG. 1. The removal of chains of oxygen atoms (colored in 
black) from the ABO3 perovskite structure (a) results in the 
vacancy ordered ABO2.5 brownmillerite structure (b), which 
consists of alternating layers of BO 4 tetrahedra and BOq 
octahedral along the b direction. The undistorted rows of 
BO4 tetrahedra can then ‘twist’ to create left- or right-handed 
chains (colored blue and red, respectively). 

eluding Ga2Al205 and Ga2Fe205).^^’^^ Each of the three 
low-symmetry bulk hettotypes display two-dimensional 
sheets of corner-connected octahedra, which rotate out- 
of-phase along the ac direction. This rotation transforms 
like the irreducible representation (irrep) Ff of the Imma 
phase. 

Along the 5 axis, however, the allowed displacements 
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FIG. 2. The hypothetical high symmetry brownmillerite struc¬ 
ture (a) is defined as having no octahedral or tetrahedral 
rotations, and has the Imma space group. Relative ordering 
of tetrahedral chains results in three different low-symmetry 
structures. If all chains are of the same handedness, the struc¬ 
ture is polar I2hm (b); alternation of left- and right-handed 
chains within each layer results in centric Pbcm (c), while 
alternation between each layer gives centric Pnma (d). The 
^-site cations are omitted from the low-symmetry structures 
for clarity. 


of the oxygen atoms are constrained by the ‘handedness’ 
and ordering of the tetrahedral chains, which ultimately 
control the final symmetry of the brownmillerite struc¬ 
ture (see the discussion in Section III). The tetrahedral 
chains can cooperatively rotate in a variety of ways, each 
described by a different irrep of Imma. When the tetra- 
hedra rotate into either all left- or all right-handed chains 
the structure displays the polar space group I2hm ow¬ 
ing to the irrep (Figure 2 b). If there is a racemic 
mixture of both types of chains, the structure becomes 
centrosymmetric, with different relative orderings gener¬ 
ating different symmetries. Alternating chains of different 
handedness within each tetrahedral layer is described by 
the A4 irrep and yields the centrosymmetric Pbcm struc¬ 
ture (Figure 2 c), while alternation between layers (given 
by X^) gives the centric Pnma structure (Figure 2 d). 
Because the left- and right-handed chains are related by 
symmetry and differ only by small atomic displacements, 
the formation energies for the different polymorphs are 
nearly degenerate and should form with equal probability. 
However, each of the aforementioned ordering types is 
seen experimentally in various members of the brownmil¬ 
lerite family, and the driving force behind the preferred 
type in different chemistries is not completely understood. 

Another degree of freedom appears for the brownmil¬ 
lerite structures when when they are grown as a thin 




FIG. 3. When brownmillerite structures are placed under epi¬ 
taxial strain, the oxygen deficient layers can order (a) parallel 
or (b) perpendicular to the substrate, with the pseudocubic 
orientations shown in (c) and (d), respectively. 


film, e.g.^ via molecular beam epitaxy or pulsed laser 
deposition, owing to the constraints imposed by epitaxial 
strain. In this thin film case, the oxygen-deficient lay¬ 
ers can order parallel or perpendicular to the substrate 
(this is shown in Figure 3 a and 3 b, respectively, with 
the pseudocubic orientation shown in Figure 3 c and 3 d). 
Although different strain states will stabilize one orienta¬ 
tion over the other, it is not always clear which will be 
preferred and why. In (La,Sr)Co205, for example, tensile 
strain stabilizes perpendicular ordering while compres¬ 
sive strain stabilizes parallel ordering;however, the 
opposite effect is observed in strained Ca2Fe205.^^ 

In this work, we seek to disentangle the structural 
and energetic effects operative in brownmillerite ox¬ 
ides through an investigation of bulk and epitaxially 
strained Sr2Fe205 and Ca2Fe205 using first-principles 
density functional theory calculations. The Sr2Fe205 
and Ca2Fe205 members are experimentally known to be 
stable in the ordered brownmillerite phase up to high tem¬ 
peratures (approximately 1200 K for Sr2Fe205 and 920 
K for Ca2Fe205);^^’^^ at ambient conditions, Sr2Fe205 
displays the Pbcm ordering,while Ca2Fe205 exhibits 
the Pnma structure.Additionally, they both contain 
only Fe^+ cations and display G-type antiferromagnetic 
ordering with high Neel temperatures (700 and 720 K, 
respectively).^^’^^’^^ The indirect band gap of both com¬ 
pounds is approximately 2.0 eV, owing to a F- to X-point 
transition. 

We find that, in agreement with experimental results of 
the bulk phases, the Pbcm and Pnma phase are the lowest 
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FIG. 4. The three characteristic angles of the A 2 B 2 O 5 brown- 
millerite structure (a). Rotations of the tetrahedral chains, 
of the octahedral network, and between the tetrahedral and 
octahedral layers are given by 4 >t, ^o, and $ot, respectively. 
These three angles are also sufficient to describe the case when 
the tetrahedral chains order parallel to the substrate under 
biaxial strain. In the perpendicular orientation (b), however, 
the changes in the out-of-plane lattice parameter owing to 
the strain state splits the octahedral and tetrahedral rotation 
angles. 


wave (PAW) potentials with the PBEsol functional, 
with valence electron configurations of for Ca, 

for Sr, for Fe, and 2 s^2p^ for 

A plane-wave cutoff of 500 eV and a 7 x 5 x 7 Monkhorst- 
Pack mesh was used during the structural relaxations.^^ 
We applied a Hubbard U correction of 5 eV using the 
Dudarev formalism to treat the correlated Fe 3 d states; we 
also enforced a G-type anti-ferromagnetic collinear spin 
ordering on the Fe atoms.Symmetry-adapted mode 
decompositions were performed using the ISODISTORT 
tool, part of the ISOTROPY software suite.Atomic 
structures were visualized using VESTA. 

To simulate the application of epitaxial strain by growth 
on a cubic [001] terminated perovskite substrate with a 
square surface net, we fix the a and c lattice parameters 
to be equal and allow the out-of-plane b axis and ions to 
fully relax (adopted from the approach of Ref. 32 ). The 
/ 26 m, Phcm^ and Pnma phases of Ca2Fe205 and the 
I 2 hm and Phcm phases of Sr2Fe205, in both the parallel 
and perpendicular orientation were strained from - 3 % to 
3 % in increments of 1 %. Because of the difference in size 
and equilibrium volume of the Ca2Fe2 05 and Sr2Fe205 
unit cells, application of the same percent strain results 
in different pseudocubic lattice parameters (a^c) for the 
two structures; for this reason, we report the strain in 
terms of ape rather than in terms of percentage. 


energy (and thus equilibrium) structures for Sr2Fe205 
and Ca2Fe205, respectively. Furthermore, we find that 
the parallel (or perpendicular) ordering is preferred under 
tensile (or compressive) strain; additionally, this switch 
occurs due to which ordering maximizes the intralayer 
tetrahedral chain distance under a specific strain state, 
analogously to the bulk phases. Finally, we report that the 
band gap of these materials is heavily influenced by epitax¬ 
ial strain owing to strain-induced changes to the intralayer 
bond angles. Although strain is well known to influence 
the band gap via control of the octahedral rotations in 
ABOs perovskites, the alternating octahedra-tetrahedra 
layers in the brownmillerites allow for a more complex 
coupling between the two. We show that, similar to the 
ABOs perovskites, an alteration of the B-O-B bond an¬ 
gles causes this response; unlike perovskites, however, it 
is the angle between the tetrahedrally and octahedrally 
coordinated irons which controls the band gap. The band 
gap of the compounds in the parallel orientation increases 
under increasing tensile strain, but decreases in the per¬ 
pendicular orientation, a discrepancy that arises from 
how the response of the out-of-plane lattice parameter 
to strain affects this bond angle differently in the two 
orientations. 


III. RESULTS AND DISCUSSION 
A. Polyhedral Rotation Definitions 

To characterize the brownmillerite structures, we con¬ 
sider three types of rotation angles: the angle between 
(i) the Fe04 tetrahedra (©t), (h) the FeOe octahedra 
(©o), and (hi) the tetrahedra and octahedra (©ot)- We 
report these angles as ©x = (180 — ^x)/2, where the 
definition of (Y = T, O, or OT) is shown in Figure 4 ; 
in this case, a larger angle indicates a larger distortion of 
the relative bond away from 180 °. In the bulk case, and 
the case of parallel orientation of vacancies under strain, 
these three angles are sufficient to describe the structure 
(Figure 4 a). In the perpendicular orientation, however, 
the applied strain and response of the out-of-plane lattice 
parameter now affect the tetrahedra and octahedra in 
different ways. To capture this, we divide and 
into two unique angles defined by in-plane (i.p.) and 
out-of-plane (o.o.p.) components (Figure 4 b). 

B. Bulk Phases 


II. COMPUTATIONAL METHODS 

All investigations were performed using density func¬ 
tional theory as implemented in the Vienna ab-initio Simu¬ 
lation Package (VASP).^^“^^ We used projector augmented- 


We first investigated the bulk phases of Sr2Fe205 and 
Ca2Fe205 with the three tetrahedral chain orderings 
shown in Figure 2. The main results are summarized in 
Table I. In the case of Sr2Fe205, we found that centrosym- 
metric Pbcm is the lowest energy structure, followed by 
/ 26 m, with Pnma being the highest. In Ca2Fe205, we 
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TABLE 1. The energetics and structure of bulk Sr 2 Fe 205 and Ca 2 Fe 205 . The tolerance factor (r) is dehned by Equation 2 
in the text. The energy difference between the different tetrahedral chain ordering structures is given as AF; each energy is 
given as the difference in meV between that phase and the lowest energy phase for each compound normalized to the number 
of formula units. The a, b, and c lattice parameters are given in Angstroms. The rotations of the tetrahedra, octahedra, and 
between the octahedra and tetrahedra are given in degrees as ©t, ©o, and ©ot, respectively; each is dehned in Figure 4. 
The average intralayer separation of tetrahedral chains is R (dehned in Figure 5), and the deviation in the bond lengths of an 
octahedra is given by A. The band gap (Eg) of each structure is given in eV. 






Sr 2 Fe 205 (Pt 

= 3.8 D, r= 

=0.976) 




Symmetry 

AE (meV/f.u.) 

a (A) 

b (A) 

c (A) 

©T D 

©o D 

0 

0 

S 

o 

R (A) 

A (xlO-'*) 

E„ (eV) 

I2bm 

12.6 

5.501 

15.402 

5.659 

24.48 

3.70 

15.41 

5.0977 

18.79 

2.19 

Pbcm 

0 

5.503 

15.407 

11.311 

24.54 

4.43 

15.39 

5.0988 

19.07 

2.15 

Pnma 

22.7 

5.499 

15.413 

5.659 

24.49 

3.64 

15.31 

5.0971 

19.03 

2.09 





Ca 2 Fe 205 (Pt 

= 1.7 D, r= 

=0.923) 




Symmetry 

AE (meV/f.u.) 

a (A) 

b (A) 

c (A) 

©T D 

©o D 

0 

0 

S 

o 

R (A) 

A (xlO-'*) 

E„ (eV) 

726m 

16.7 

5.381 

14.617 

5.579 

27.92 

6.538 

20.34 

5.0173 

11.16 

2.16 

Pbcm 

23.3 

5.386 

14.627 

11.155 

28.16 

4.909 

19.90 

5.0476 

10.85 

2.13 

Pnma 

0 

5.397 

14.632 

5.561 

27.23 

8.032 

20.45 

4.9941 

10.73 

2.07 


found the Pnma structure to be lowest in energy, followed 
by I2hm and Phcm. Each polymorph is separated from 
the other two by a small amount of energy (AE, Table I), 
indicating that the formation of a right- or left-handed 
tetrahedral chain may be equally probable. Furthermore, 
the lattice parameters, rotation angles, and band gaps 
vary only slightly between orderings. Between the two 
chemistries, Ca 2 Fe 205 has a smaller unit cell and larger 
rotation angles than Sr 2 Fe 205 , which is due to the smaller 
size of the Ca^+ cation. 

Experimentally, the structure of Sr 2 Fe 205 has been 
highly contested. Initial structure refinements on powder 
samples ambiguously supported assignment of both com¬ 
pletely disordered tetrahedral chains (space group Imma) 
or pure left- or right-handed ordering (space group /26m); 
this material was thus theorized to display /26m symme¬ 
try locally, but with random ordering taking place over 
longer length scales. This was challenged, however, by 
transmission electron microscopy results showing clear in¬ 
tralayer alternation of tetrahedral chains.Additionally, 
more recent neutron diffraction experiments on single 
crystals grown by the floating zone method have indi¬ 
cated that Pbcm is indeed the preferred structure type 
for Sr 2 Fe 205 .^^’^^ Furthermore, all of these recent experi¬ 
ments are supported by the aforementioned results of our 
ab initio calculations; the small energy difference between 
the Pbcm and /26m structures, however, indicate that 
intergrowths of the phases is not out of the question. In 
contrast to the ambiguity of Sr 2 Fe 205 , the structure of 
Ca 2 Fe 205 is well known to form in the Pnma phase, 
again in agreement with our theoretical results. 

What factors lead to and stabilize the preferred ground 
state in different brownmillerite compounds? In 2005, 
Abakumov et al. put forth the idea that the twisting of 
tetrahedral chains away from the undistorted 180° orien¬ 
tation creates local dipole moments, with larger rotations 
producing larger dipoles.Hadermann and Abakumov 
et al. further suggested that the distance (z.e., the length 


of the 6 axis) between the tetrahedral layers is also an 
important factor to consider.Parsons et al. then rec¬ 
ognized that each tetrahedral ordering scheme distorts 
the octahedra in different ways; the fact that the octa¬ 
hedra are not connected out-of-plane (along 6) causes 
the apical oxygen atoms to displace more than the equa¬ 
torial ones, which creates a ‘shearing’ effect from this 
non-rigid rotation. Generally, the /26m phase causes the 
least octahedral (elastic) distortion, followed by Pnma, 
with Pbcm causing the most. They then rationalized 
these arguments into a “structure map” relating the ob¬ 
served phase of different brownmillerites to these fac¬ 
tors, followed by a classification of several known com¬ 
pounds into this scheme.Although these observations 
have been key in building an understanding of structural 
trends in brownmillerite oxides, and are corroborated by 
some recently synthesized brownmillerite phases (such 
as Ca 2 Cr 2 © 5 ),^^discrepancies in this structure map show 
there are additional effects which should be considered. 
Ca 2 FeCo 05 and Ca 2 Co 205 , for example, both display 
the Pbcm rather than the structure-map predicted Pnma 
structure. 

We now seek to explain the stability of these materials’ 
preferred ground state phase in terms of local structure. 
The tetrahedral chain ordering configuration preferred 
by the ground state can be summarized as a complex 
competition between two primary energetic factors: (z) 
separation of the tetrahedral chains (z.e., minimization of 
electrostatic repulsion) and (ii) distortion of the nominally 
regular BOq octahedra (minimization of elastic strain 
energy). The three different ordering schemes shown in 
Figure 2 better maximize either factor (z) or (zz), at the 
expense of the other. Following the approach of Zhang 
et al.^^^ we can assign a magnitude to the local dipole 
generated by the tetrahedral rotations (Pt, Table I, given 
in units of Debye). In compounds with large Pt, factor (z) 
becomes the quantity of interest to maximize; rather than 
considering the separation along 6, however, we quantify 
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FIG. 5. The intralayer chain separation in the brownmillerite 
structure. The average separation distance is given by R, 
which is dehned as the average of the two distances between 
Fe atoms of different chains {di and 0 ^ 2 ). 


this by defining the average m^ralayer separation of the 
tetrahedral chains as R (shown in Figure 5). Competing 
with this is factor (zz), z.e., the regularity in the octahedra. 
We use the averaged sum-of-squares difference between 
the measured bond lengths (d^) and the average bond 
length (davg) in the octahedra to quantify this: 



For small values of Pt, factor (zz) becomes the critical 
quantity to optimize, and minimizing A becomes more 
important than maximizing R. 

Interestingly, we found that although Sr 2 Fe 2 05 has 
smaller rotations, it has a larger chain dipole (Pt = 3.8 
D) than Ca 2 Fe 2 05 (Pt = 1.7 D); generally, one would 
expect smaller rotations to generate a smaller dipole. This 
is due to the fact that the electric polarization generated 
by the apical oxygen atoms tends to cancel that from 
the equatorial oxygen atoms in the tetrahedra; however, 
the smaller rotations of the octahedra in Sr 2 Fe 205 means 
this cancellation happens to a lesser degree, increasing 
Pt- Again, at large values of Pt, the need to separate 
the dipoles, factor (z), overcomes the need to have regular 
octahedra, factor (zz). Sr 2 Fe 205 therefore exhibits the 
Pbcm structure, which allows for the best average in¬ 
plane separation (P, Table I) of the large dipoles. In 
Ca 2 Fe 205 , the larger octahedral rotations distort the 
octahedra more in the presence of tetrahedral chains; this, 
in combination with the small chain dipole, therefore gives 
more importance to the phase which keeps the octahedra 
most regular (A, Table I). In this case, that is the Pnma 
structure. 

Another approach to understanding these energetic 
competitions is to consider the the Goldschmidt tolerance 
factor.This quantity is given by: 


rAprp 

V2{rB + ro) ’ 


( 2 ) 


where tb^ and rp is the radius of the A-site atom, 
P-site atom, and oxide anion, respectively.^^ An ideal 
undistorted (cubic or tetragonal) structure will have r=l; 


the further the tolerance factor deviates from 1, the more 
the structure is likely to distort, and the larger the oc¬ 
tahedral rotations will be. In Ca 2 Fe 2 05 (r=0.923), the 
Pnma phase least distorts the octahedra (as refiected in 
A in Table I), and is therefore the most stable. As r 
increases, the need to maximize the separation between 
dipoles overcomes the steric packing driving force for 
octahedral rotations and hence they decrease; because 
Pbcm best maximizes this distance (P in Table I), it is 
the preferred phase for the larger Sr 2 Fe 205 (t=0.976). 

Finally, despite the structural differences between the 
bulk Sr 2 Fe 205 and Ca 2 Fe 205 phases, the DFT-PBEsol 
(indirect F-X) band gap remains nominally the same 
(^ 2.10 eV) across changes in chemistry and tetrahedral 
ordering. This is due to the fact that O 2 p states make 
up the top of the valence band, while the bottom of the 
conduction band is made up of Fe 3d states, making the 
band gap largely independent of A-site chemistry. A more 
detailed analysis of the changes in the electronic gap is 
given for the thin film cases below. 


C. Strained Phases 

We next investigated the crystal and electronic struc¬ 
ture of Sr 2 Fe 205 and Ca 2 Fe 205 under epitaxial strain. 
As mentioned previously, we must now consider the orien¬ 
tation of the vacancies as another degree of freedom (z.e., 
whether the vacancy layers are parallel or perpendicular to 
the substrate. Figure 3), in addition to tetrahedral chain 
rotation and ordering. First, we find that when Sr 2 Fe 205 
is constrained to a substrate, it retains the Pbcm ground 
state in both parallel (filled symbols) and perpendicular 
(empty symbols) substrate-vacancy orientations (Figure 
6a, top), a fact that can be explained using the same 
arguments presented for the bulk phase. However, we 
also find that the perpendicular configuration of oxygen- 
deficient layers is stabilized over the parallel arrangement 
under small amounts of compressive strain. This is due 
to the fact that increasing compressive strain decreases 
the intralayer separation between tetrahedral chains if 
the vacancies are ordered parallel to the substrate (ie., 
becomes more energetically unfavorable), but increases 
R if the vacancies are perpendicular. The strain value at 
which the perpendicular ordering becomes more stable 
exactly coincides with the values at which the intralayer 
separation becomes greater (Figure 6a, middle) and the 
octahedral distortions become smaller (Figure 6a, bottom) 
than that of the parallel ordering. 

One important detail to note is that Sr 2 Fe 2 05 has a 
relatively large pseudo-cubic lattice parameter {ape); only 
with some of the largest commercially available substrates 
can it be placed under tensile strain and have the parallel 
orientation be stabilized. Many experimental observa¬ 
tions agree with our prediction of the energetic stability 
of vacancy orientation. Growth of Sr 2 Fe 205 near the 
transition point, such as on SrTiOs {ape = 3.91 A) shows 
a competition between the two orientations,whereas 
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FIG. 6. Energy of the different tetrahedral chain ordered Sr 2 Fe 205 (left) and Ca 2 Fe 205 (right) structures as a function of 
epitaxial strain (top panels). In both cases, the parallel vacancy ordering (filled symbols) is stabilized under tensile strain, while 
perpendicular (empty symbols) is stabilized under compressive. This change in stabilization occurs at the point where either the 
parallel or perpendicular phase maximizes the average intralayer tetrahedral chain separation (R, middle panels) and minimizes 
the octahedral distortion effect (A, bottom panels). 


growth on larger substrates stabilizes the parallel orienta¬ 
tion, such as on KTaOs = 3.99 A)."^^ 

When Ca 2 Fe 2 05 is placed under epitaxial strain, the 
Pnma phase remains lowest in energy, while the I2hm 
and Phcm structures become much closer in energy and 
strongly compete. As with Sr 2 Fe 205 , the tetrahedral lay¬ 
ers in Ca 2 Fe 205 switch to a perpendicular orientation un¬ 
der compressive strain, and parallel under tensile (Figure 
6b, top). Once again, this occurs due to the perpendicular 
orientation maximizing the distance R between intralayer 
tetrahedral chains and minimizing octahedral distortions, 
A, under compression (Figure 6b, middle and bottom). 
Although the Pnma phase minimizes octahedral shearing 
as in the bulk phase, the very small differences in A be¬ 
tween the strained phases means that it cannot be the only 
controlling feature; instead there is a complex interplay 
between the various structural descriptors which produces 
the observed ground state. The fact that Ca 2 Fe 205 has 
a smaller pseudo-cubic lattice parameter than Sr 2 Fe 205 
means that the parallel phase is experimentally accessible 
with much more modest strain. Experimental results on 
thin films of Ca 2 Fe 205 show that the vacancies order 
perpendicularly when grown on LaSrA 104 (LSAO, ape = 
3.75 A) and LaAlOs (a^c = 3.79 A), but order parallel on 
LSAT {ape = 3.87 A) and SrTiOs {ape = 3.91 A), again 
in agreement with these theoretical results.Interest¬ 
ingly, the stabilization of the perpendicular phase under 
compressive strain and parallel phase under tensile strain 
does not hold across all brownmillerite oxide chemistries; 
in the cobaltates, for example, the trend is reversed. 

We next look at the structural evolution of these com¬ 
pounds under strain. For the structures with the parallel 
orientation, there are only the three bond angles defined 
previously to consider (Figure 4). Under increasing tensile 


strain, the Fe 04 tetrahedral chains become less distorted 
(decreasing Fe-O-Fe bond angle), while the octahedral 
network and connections between the tetrahedral and 
octahedral layers, becomes more distorted (increasing Fe- 
O-Fe bond angle). It should come as no surprise that 
this trend holds across all three types of ordered phases, 
as they are all very structural similar with the exception 
of the relative orientation of tetrahedral chains (Figure 
7). Additionally, the chemistry of the A-site, although 
influencing the magnitude of the tilts, does not affect this 
trend. 

When the oxygen-deficient layers reorient to become 
perpendicular to the substrate, there are now two more 
angles to consider. In the previous geometries, only one 
angle of the tetrahedra and octahedra had to be considered 
owing to the fact that the equatorial oxygen atoms of 
each were in the plane of epitaxial strain. Now these same 
atoms are oriented such that they are affected by both the 
plane of epitaxial strain and the strain-induced changes to 
the out-of-plane lattice parameter. In order to get a full 
understanding of the structural distortions in this case, 
we separate the tetrahedral and octahedral angles into in¬ 
plane and out-of-plane components. In this perpendicular 
orientation, increasing tensile strain increases the out-of¬ 
plane component of the octahedral and tetrahedral angles, 
while decreasing the in-plane component, as well as the 
angle between the layers. 

Finally, the band gap of both materials is strongly 
influenced by strain, ranging from 1.8 eV to 2.6 eV for 
Sr 2 Fe 205 (Figure 8a) and 1.6eV to 2.5eV for Ca 2 Fe 205 
(Figure 8b). In both compounds, increasing tensile strain 
results in an increase of the band gap for the parallel ori¬ 
entation of vacancies, but interestingly, a decrease in the 
perpendicular orientation. This is due to how an increase 
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FIG. 7. The effect of strain on the tilt angles (as defined in Figure 4) of the different Sr 2 Fe 205 (left) and Ca 2 Fe 205 (right) 
structures with both parallel and perpendicular vacancy orientation. 


in the lattice parameters influences the connectivity of 
the tetrahedrally and octahedral coordinated iron atoms 
(z.e., ©ot)- For the films with the parallel vacancy orien¬ 
tation, ©OT deviates further away from 180 ° as tensile 
strain increases; this results in decreased overlap of the Fe 
d-orbitals and O p-orbitals, giving a higher band gap. Fig¬ 
ures 8c and 8d show, as an example, the difference in the 
band gap of parallel-oriented Sr2Fe205 at 2% ( 2.43 eV) 
and -2% ( 2.00 eV). The opposite effect occurs in the per¬ 
pendicular case, where a decrease in ©ot provides better 
overlap and thus a smaller band gap through increased 
bandwidth. 


D. Comparison to Perovskite Oxides 

The fact that BOq octahedra in perovskite oxides form 
a flexible corner-connected network (Figure la) allows 
them to easily rotate in space about the different crys¬ 
tallographic axes. The size of these rotations directly 
affect the magnitude of the B-O-B bond angles, which 
in turn impacts many electronic and magnetic proper¬ 
ties. Although there are 15 distinct ways in which the 
octahedra can cooperatively rotate while retaining con¬ 
nectivity (as identified by Glazer),"^^ the vast majority 
exhibit either an orthorhombic or rhombohedral tilt pat¬ 
tern (given by a~a~c^ or a~a~a~ in Glazer notation, 
respectively). A high degree of control over the electronic 
structure can be achieved by using epitaxial strain (or 
chemical substation, as captured by r) to control these 
rotations, owing primarily to the strong coupling between 
the lattice and electronic degrees of freedom in perovskites 
(Figure 9 a);^^“'^^ as mentioned previously, buckling of the 
B-O-B bond away from a linear 180 ° configuration (z.e.. 


increasing the magnitude of the octahedral rotations) 
decreases the overlap between the O p and metal B d 
orbitals, thereby increasing the band gap in insulating 
compounds or inducing bandwidth-driven metal-insulator 
transitions. 

As we have shown, the distinct alternating tetrahe¬ 
dral and octahedral layers in brownmillerites allows for 
many more structural degrees of freedom than perovskites. 
Application of the same forces (chemical pressure and 
epitaxial strain) now affects the interlayer separation of 
tetrahedra (ie., b lattice parameter), as well as the ro¬ 
tations (connectivity) of the different polyhedra (given 
by ©o, ©T, and ©ot); each of these in turn influences 
the octahedral distortions (A), the intralayer dipole sep¬ 
aration (R ), and the magnitude of the local tetrahedral 
dipoles (Pt)- The balance of these different factors, which 
is shown schematically in Figure 9 b, is then what governs 
the equilibrium structure—the tetrahedral chain order¬ 
ing and, for the case of thin films, the orientation of the 
vacancies relative to the substrate. 

A feature implicitly indicated in Figure 9 is that the 
electronic properties of these oxides are also a structural 
consequence; hence, control over the structural descrip¬ 
tors makes it possible to control the electronic response. 
Indeed, the band gaps of the brownmillerites are highly 
sensitive to rotations of the polyhedra; like the perovskites, 
the gap opens as the rotations increase owing to a change 
in orbital overlap. Unlike perovskites, however, it is not 
solely a single B-O-B bond angle controlling the elec¬ 
tronic structure, which is due to the breaking of the 
structural topology by the oxygen-vacant (tetrahedral) 
layers. For the brownmillerite oxides, it is the out-of-plane 
angle between the PO4 tetrahedra and BOq octahedra, 
given by ©ot (or, alternatively, Ptet“0“Poct), that con- 
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FIG. 8. The band gap of (a) Sr 2 Fe 205 and (b) Ca 2 Fe 205 are strongly influenced by epitaxial strain; the electronic gap for 
thin films structures with vacancies ordered parallel to the substrate (hlled symbols) increases as tensile strain increases, but 
decreases in the perpendicular orientation (empty symbols). This electronic structure response occurs owing to the manner in 
which strain affects the angle between the tetrahedral and octahedral layers. The Pbcm structure of Sr 2 Fe 205 with vacancies 
parallel to the substrate, for example, has a ~ 0.3 eV larger band gap under (c) 2% tensile strain when compared to (d) 2% 
compressive strain as seen by the change in the atom-resolved densities of states. 
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FIG. 9. The effect of different stress stimuli (ionic size or chem¬ 
ical pressure [captured by r] and epitaxy) on various structure 
descriptors, which combine to produce the equilibrium crystal 
structure in (a) perovskite and (b) brownmillerite oxides. 


trols the gap and band edge character. Furthermore, 
strain affects the band gap of parallel or perpendicularly 
oriented brownmillerites in completely opposite (asym¬ 
metric) ways owing to the response of ©ot to the biaxial 
strain state. Although strain has been used to control 
functional properties in perovskites, it appears to induce 
even more interesting responses in brownmillerites owing 
to the additional structural degrees of freedom present. 


IV. CONCLUSION 

Using first principles density functional theory calcu¬ 
lations, we investigated the brownmillerite compounds 
Sr 2 Fe 205 and Ca 2 Fe 205 . These anion deficient materials 
offer many more degrees of freedom than the related per¬ 
ovskite structures, increasing the structural complexity 
greatly. The symmetry of these structures is determined 
by the relative ordering of left- and right-handed tetrahe¬ 
dral chains; which ordering scheme preferred, however, is 
determined by the interplay of several structural factors. 
Although previous studies have attempted to characterize 
the various members of this family, here we clearly demon¬ 
strate how the tolerance factors, polarizations generated 
by rotations of Fe 04 tetrahedra, distortions of FeOe octa- 
hedra, and intralayer chain separation interact to produce 
the observed equilibrium phases. 

Furthermore, we show the crystal and electronic struc¬ 
ture of these phases are highly tunable using epitaxial 
strain. We find that under compressive strain, the pre¬ 
ferred ordering of vacancies (tetrahedral Fe 04 chains) 
changes from parallel to perpendicular when the separa¬ 
tion between the change becomes better maximized. The 
rotation angles are also strongly influenced, resulting in 
a increase (or decrease) in the band gap for the paral¬ 
lel (or perpendicular) orientation. While this study is 
by no means comprehensive in scope, we hope that the 
framework presented here encourages the consideration of 
additional structural metrics (especially intralayer sepa¬ 
ration) in future studies and helps provide greater insight 
into this complex family of materials. We can now harness 
this knowledge and understanding of the atomistic mech¬ 
anisms behind the coupling of the crystal and electronic 
structure in brownmillerites to design electronic function 
and realize new materials platforms for thin film devices. 
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